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There are several ways to create the vorticity-free solitary waves - rarefaction pulses - in conden- 
sates: by the process of strongly nonequilibrium condensate formation in a weakly interacting Bose 
gas, by creating local depletion of the condensate density by a laser beam, and by moving a small 
object with supercritical velocities. Perturbations created by such waves colliding with vortices are 
studied in the context of the Gross-Pitaevskii model. We find that the effect of the interactions 
consists of two competing mechanisms: the creation of vortex line as rarefaction waves acquire 
circulation in a vicinity of a vortex core and the loss of the vortex line to sound due to Kelvin 
waves that are generated on vortex lines by rarefaction pulses. When a vortex ring collides with 
a rarefaction wave, the ring either stabilises to a smaller ring after emitting sound through Kelvin 
wave radiation or the entire energy of the vortex ring is lost to sound if the radius of the ring is 
of the order of the healing length. We show that during the time evolution of a tangle of vortices, 
the interactions with rarefaction pulses provide an important dissipation mechanism enhancing the 
decay of superfluid turbulence. 

PACS numbers: 03.75.Lm, 05.45.-a, 67.40.Vs, 67.57.De 



I. INTRODUCTION 

Low temperature superfluids and recently discovered 
dilute Bose-Einstein condensates are often modelled by 
the Gross-Pitaevskii equation (GP), which is a nonlin- 
ear Schroedinger equation on the one-particle wave func- 
tion, ip. This equation provides a simple framework for 
study many fundamental hydro dynamical properties of 
condensates. In particular, vortex-sound interactions in 
condensates are receiving increasing attention over last 
couple of years fl 0, H H H, • It has been suggested 
that the emission of sound by vortex reconnections and 
vortex motion is the only active dissipation mechanism 
responsible for the decay of superfluid turbulence. The 
decay of superfluid turbulence via Kelvin wave radiation 
and vortex reconnections was studied in the framework 
of the GP equation . In this study the collision of two 
small vortex rings were analysed and the loss of the vor- 
tex line due to reconnections and Kelvin wave radiation 
was numerically evaluated. 

The goal of this paper is to consider the effect of the 
vorticity-free solitary waves that together with vortices 
are created in Bose-Einstein condensates either during 
the self-evolution of a Bose gas from a strongly nonequi- 
librium initial state or evolved from local density deple- 
tions of a condensate. We show that the interactions with 
such waves trigger the loss of the vortex line via Kelvin 
wave radiation and that these interactions enhance the 
dissipation of the vortex tangle. 

We write the GP equation in dimensionless form as 
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in dimensionless variables such that the unit of length 
corresponds to the healing length £, the speed of sound is 



c = 1/V% and the density at infinity is = l^ool 2 = 1- 
The numerical integration was performed using a finite 
difference scheme. The faces of the computational box 
were open to allow sound waves to escape; this is achieved 
numerically by applying the Raymond-Kuo technique Q • 
In the turbulence simulations the periodic boundary con- 
ditions were used to conserve the total number of parti- 
cles in the system. 

In a seminal paper, Jones and Roberts numerically 
integrated the GP equation (QJ and determined the entire 
sequence of solitary wave solutions of the GP equation, 
such as vortex rings, vortex pairs, and finite amplitude 
sound waves named rarefaction pulses. They showed the 
location of the sequence on the momentum, p, energy, 
£ , plane. In three dimensions they found two branches 
meeting at a cusp where p and £ assume their minimum 
values, p m and £ m . As p — > oo on each branch, £ — > oo. 
On the lower branch the solutions are asymptotic to large 
vortex rings. 

As £ and p decrease from infinity along the lower 
branch, the solutions begin to lose their similarity to large 
vortex rings. Eventually, for a momentum po slightly 
greater than p m , they lose their vorticity loses its 
zero), and thereafter the solitary solutions may better be 
described as 'rarefaction waves'. The upper branch con- 
sists entirely of these and, as p —> oo on this branch, the 
solutions asymptotically approach the rational soliton so- 
lution of the Kadomtsev-Petviashvili Type I equation. 

The Jones-Roberts (JR) solitons are the only known 
disturbances that propagate with a constant velocity. 
Notice, however, that there many other waves that 
change their velocity and shape during their motion. For 
instance, a strong perturbation of a rarefaction pulse on 
the upper branch of the dispersion curve causes it to col- 
lapse onto the lower branch and become a vortex ring. 
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During its transition the wave looses its energy and mo- 
mentum and the minimum of its density decreases grad- 
ually before it reaches zero. In the next stage of the evo- 
lution the wave becomes a non-axisymmetric ring which 
radius increases until it reaches an axisymmetric solitary 
state on the lower branch of JR dispersion curve. All 
these intermediate time snapshots of the wavefunction of 
the condensate have higher energy than any solitary wave 
on the lower branch below the final axisymmetric vor- 
tex ring state although they may have the same minima 
of the density in the transition. In what follows we re- 
serve the term "rarefaction pulse" or "rarefaction wave" 
to describe a finite amplitude sound wave which is a JR 
solitary wave. The scenario just described illustrates a 
typical mechanism in which the vortex line length may 
be created and increased as the result of perturbations 
of rarefaction pulses. In Q we have also considered a 
creation of vortex rings as a result of energy and momen- 
tum transfer between two interacting rarefaction pulses. 
But given a complex tangle of interacting vortices we ex- 
pect that the loss of the vortex line length will dominate 
its creation to account for the experimentally observed 
decay of superfluid turbulence. 

In |?j we have developed an algorithm for finding ap- 
proximations to the JR solitary wave solutions with the 
correct asymptotic behaviour at infinity. An axisym- 
metric solitary wave moving with the velocity U along 
the x— axis is accurately approximated by ip(x, s) = 
1 + u(x, s) + iv(x, s) where 
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where a^- , bij , Cij and the dipole moment m can be de- 
termined from the series expansion and are functions 
of U. Also in (0), s 2 = y 2 + z 2 , x = x - Ut, and 
f(U) = 1 — 2U 2 . Notice, that (0) represents a vortex 
ring as well as a vorticity-free rarefaction pulse depend- 
ing on U. If — 2 < aoo < — 1, then J2J) represents a vortex 
ring as the power series expansion around zero shows; if 
— 1 < aoo < 0, then (0) gives an approximation of a rar- 
efaction pulse, and aoo = — 1 is a borderline case, such 
that the solitary wave has a single zero of the wavefunc- 
tion and, therefore, can be called a point defect. These 
approximations can be used as initial conditions in the 
numerical simulation that study interactions, when the 
initial state is prepared by multiplying the wavefunctions 
of the distant individual solitary waves. Without an ac- 
curate starting point in numerical calculations it would 
be impossible to separate clearly the effect of interac- 
tions from the evolution of each solitary wave by itself as 
it settled down from a poor initial guess. 



II. CREATION OF RAREFACTION PULSES IN 
CONDENSATES 

Here we study the excitations created by the collisions 
between rarefaction pulses and vortices and the effect 
these collision have on the evolution of vortex line length 
in a regime of superfluid turbulence with a tangle of vor- 
tices. We start by describing three typical scenario in 
which a large amount of rarefaction pulses is created, so 
that their collisions with vortices do become significant. 
The process of strongly nonequilibrium BEC formation in 
a macroscopically large uniform weakly interacting Bose 
gas was elucidated in using numerical integration of 
the GP equation. As the system evolves from weakly 
turbulent state to state of strong turbulence, the phases 
of the complex amplitudes of the wave field ip become 
strongly correlated and the period of their oscillations 
become comparable with the evolution times of the oc- 
cupation numbers. This signifies that the quasiconden- 
sate is formed with appearance of a well-defined tangle 
of quantized vortices and localized vorticity-free solitary 
structures such as rarefaction pulses manifesting the start 
of the final stage in the Bose gas evolution: the decay of 
superfluid turbulence. The right panel of Fig. 5 pre- 
sented in [HI shows a single vortex ring surrounded by 
many rarefaction pulses as a result of the turbulent decay 
of the initial vortex tangle, which also implies that many 
more rarefaction pulses were present at the beginning of 
the decay. 

Next we describe how the JR solitons and in partic- 
ular rarefaction pulses can be created in Bose-Einstein 
condensates. It is generally believed that to create vor- 
tices it is necessary to transmit angular momentum by 
rotationally stirring the condensate with a laser beam. 
In [m we demonstrated that the collapse of a stationary 
spherically symmetric bubble can lead to the vortex nu- 
cleation. After the condensate fills the cavity, it begins 
to expand with growing density oscillations. These 'dips' 
of the density are themselves unstable and the develop- 
ment of this instability leads to the creation of localized 
disturbances: vortex rings and rarefaction pulses. 

Similarly, the evolution of a depletion of the condensate 
by laser beam after the laser was turned off may lead to 
creation of JR solitons. Notice, that there is no need 
to add angular momentum or deplete the density of the 
condensate to zero. We demonstrate this by considering 
a depletion of the condensate, such that the wavefunction 
of the condensate just before the laser is switched off is 
given by 



V(x) = 



\ tanh 



0.01(> 2 + 0.1(j/ 2 



z 2 ) 



100) . (3) 



This density depletion roughly corresponds to the exper- 
iments in the sodium condensates, see for instance [l2j . 
with the healing length £ = 0.3/im and the Gaussian 
beam waist of more than 10£. The oblate spheroidal 
form of the depletion (J3J would be formed if such a laser 
beam was guided in the condensate in a short straight 
streak of the length of about 50£. 
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FIG. 1: (colour online) The snapshots of the contour plots of 
the density cross-section of a condensate obtained by numer- 
ically integrating the GP model (QJ for the oblate spheroidal 
density depletion Black solid lines show zeros of real 

and imaginary parts of ip, their intersection shows the posi- 
tion of topological zeros. Both low and high density regions 
are shown in darker shades to emphasise intermediate density 
regions. 
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We integrated the GP equation (QJ using Q as our 
initial field in the box of the volume V = 120 3 . The 
faces of the computational box were open to allow sound 
waves to escape. Fig. 1 shows the time snapshots of the 
density cross-sections for z = 0. Twelve vortex rings of 
various radii and four rarefaction pulses moving outward 
from the center of the depletion are clearly seen at t = 
73 snapshot. Fig. 2 gives the isoplots of the density at 
|^| 2 = 0.3. Rarefaction pulses are seen as small oblate 
spheroids. 

The nucleation of vortices in a uniform condensate has 
been linked to critical velocities of the flow [ill 13 • 
It has also been pointed out that a moving object of 
a size comparable with the healing length generates rar- 
efaction pulses rather than vortex rings when the velocity 
on the surface of this object exceeds the local speed of 
sound. Therefore, it is possible to generate the rarefac- 
tion pulses by guiding small objects with supercritical 
velocities through the condensate. 



FIG. 2: (colour online) Time snapshots of the density iso- 
plots p — 0.3 of the evolution of the initial depletion of the 
condenstate amplitude given by @. 
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III. COLLISIONS WITH A STRAIGHT LINE 
VORTEX 

In our first example we consider a collision between 
a rarefaction pulse and the straight-line vortex of unit 
winding number. The rarefaction pulse is moving with 
a constant self-induced velocity U = 0.63 and belongs 
to the lower branch of the JR dispersion curve. The 
sequence of density isosurface plots illustrating the col- 
lision for two impact angles (tt/2 and 7r/4) is shown in 
Fig. 3. Initially the center of the rarefaction pulse is 
20 (Fig. 3(a)) or 10 (Fig. 3(b)) healing lengths apart 
from the axis of the vortex. During the collision the rar- 
efaction pulse creates a distortion on the vortex line by 
exciting two Kelvin wave packets. In the close vicinity of 
the vortex line, the rarefaction pulse acquired vorticity it 
previously did not have and became a vortex ring (Fig. 
3(a) at t = 25 and t = 30 and Fig. 3(b) at t = 20), so we 
would expect that the result of the collision is similar to 
what was found in _6] for collisions of large vortex rings 
with small rings. 

The dispersion relationship for Kelvin wave in the GP 
model in the limit k — > is [l7| 

u ~ \k 2 (\n(k) + 0.003187...), (4) 

which gives a group velocity of the Kelvin wave packet 
as 

C/ s = -^21n[A]_ 1.00637^), (5) 

where A is the dominating wavelength. Our calcula- 
tions show that A ~ 8 healing lengths and U g ~ 0.2 
(U g ~ 0.3c), where c is the speed of sound) which fits (0 
quite well. This result also agrees (when different scal- 
ing of the GP equation is taken into account) with the 
wavelength of the excited Kelvin wave during the vortex 
reconnections in 

The collision with a rarefaction pulse from the upper 
branch of the JR dispersion curve excites a Kelvin wave 




FIG. 3: Sequence of density isosurfaces for p — 0.3 illustrat- 
ing a collision of the rarefaction pulse moving with the velocity 
U — 0.63 and the straight line vortex. Two impact angles are 
shown: tt/2 in (a) and 7r/4 in (b). The collision excites a pair 
of Kelvin wave packets propagating in the opposite directions 
along the vortex line. 
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packet of a similar central wavelength. Fig. 4 shows the 
contour plots of the density cross-sections of two such 
collisions for two different rarefaction pulses that belong 
to lower (left panel) and upper (right panel) branches of 
the JR dispersion curve. The wavelength of the created 
excitation is best seen by the intersections of zeros of 
real and imaginary parts of the wavefunction ip shown 
by solid and dashed lines correspondingly. Similarly, a 
collision with an offset generates a Kelvin wave packet 
of approximately the same dominating wavelength but a 
smaller amplitude. 

We can estimate the amplitude of the Kelvin wave gen- 
erated by a distant rarefaction wave by reducing the 
problem to two dimensions. We shall assume that ini- 
tially the straight line vortex is positioned at (0, 0, z) and 
the rarefaction wave is moving in the positive x— direction 
along y = —d y < with £ <C d y . Under these assump- 
tions, the wavelength of the perturbation along the vor- 
tex line is much larger than its amplitude, 77, (the maxi- 
mum displacement of the vortex center from the origin in 
xy— plane), therefore, we can assume that the x— and y— 
components of the force exerted by a rarefaction wave on 
the vortex dominate the z— component (in other words, 



FIG. 4: (colour online) The snapshots of the contour plots 
of the density cross-section of a condensate obtained by nu- 
merically integrating the GP model © illustrating the col- 
lision of the rarefaction pulses and the straight line vor- 
tex at t — 35 that are initially a distance of 20 healing 
lengths apart. Left panel shows the collision of the rarefac- 
tion pulse from the lower branch of the JR dispersion curve 
(U = 0.63,8 = 52.3, p = 72.2, m = 8.37) and the right panel 
is for the rarefaction pulse from the upper branch of the JR 
dispersion curve (U — 0.68, £ = 53.7, p = 74.1, m — 8.8). 
Black solid and dashed lines show zeros of real and imaginary 
parts of ip correspondingly, therefore, their intersection shows 
the position of topological zeros. Both low and high density 
regions are shown in darker shades to emphasise intermediate 
density regions. Only a portion of an actual computational 
box is shown. 
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the vortex may be considered as rigid). Also, if the inci- 
dent angle of the rarefaction pulse and the direction of the 
vortex axis is nonzero, the amplitude calculated below 
would have to be multiplied by the cosine of this angle. 
In xy— plane the vortex is advected by the rarefaction 
pulse, so that it instantaneous velocity coincides, to the 
leading order, with the velocity of the rarefaction pulse at 
the vortex location. The center of the rarefaction pulse is 
moving with the self-induced velocity, U = (£7, 0, 0) and 
is advected by the vortex. The motion of the vortex and 
of the rarefaction pulse is, therefore, given by the system 
of the coupled ordinary differential equations 

dx^_ y r - y v 

dt ~ (x r - x v ) 2 + (y r - y v ) 2 ' 

dy r x r - x v 

dt ~ (x r - x v ) 2 + (y r - y v ) 2 ' ^ ] 

dr v 

— = V x (x v -x r ,y v -y r ), 
dv v 

-J- = V y (x v -x r ,y v -y r ), 

where (x r (t), y r (t)) and (x v (t) , y v (t)) are the positions 
of the centres of the rarefaction pulse and the vortex, 
correspondingly, at time t. V x and V y are the velocity 
components induced by the rarefaction pulse on the vor- 
tex given by V x = dS/dx and V y = dS/dy, where S is 
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FIG. 5: (colour online) The parametric plots of the trajecto- 
ries (x v (t), y v (t)) of the vortex initially at (0, 0) and moving in 
xy— plane under the influence of the rarefaction pulse initially 
at (-20,-20), (-20,-10), (-10,-20) and (-10,10). The 
vortex travelled the distance between two adjacent points in 
unit time. 



FIG. 6: The amplitude of the maximum displacement of the 
vortex from the origin as the result of the interactions with 
the rarefaction pulse from the lower (U = 0.63) and upper 
(U — 0.69) branches of the JR dispersion curve as functions 
of the initial distance d y between the vortex and the center 
of the rarefaction pulse along the y— axis. 
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the phase of the wavefunction i/j = Rexp lS , so that 
(l+u)udv-vdu 

— — — r^— ; 2 ' \' ) 

(1 + u) z + V A 

where the real, 1 + u, and imaginary, v, parts of ip are 
approximated by (|2J). The initial conditions are given 
by (x v (0),y v (0)) = (0,0) and (x r (0) , y r (0)) = -d = 
— (d x , d y ). Fig. 5 shows the trajectories of the vortex for 
U = 0.63 and d = (10, 10), (10, 20), (20, 10) and (20, 20). 
It is clear from the figure, that the maximal amplitude 
of the displacement, 77, does not depend strongly on d x , 
but on dy. Fig. 6 plots 77 as a function of dy for two rar- 
efaction pulses that belong to the lower (U = 0.63) and 
upper (U = 0.69) branches of the JR dispersion curve. 
For large offsets, d y , (d y > 30) the maximum amplitude, 
77, decays exponentially with d y , so that 77 ~ d y - 24 for 
U = 0.63 and 77 - d^ 62 for U = 0.69. 

IV. COLLISIONS WITH A VORTEX RING 

Next we consider the collision of a rarefaction pulse 
with a vortex ring. As in the previous cases, such a col- 
lision excites Kelvin wave with release of sound energy. 
Fig. 7 illustrates the time snapshots of the collision be- 
tween the vortex ring of radius R = 5.1 moving along the 
x— axis and the rarefaction pulse moving along the line 
y = —R. The total vortex line is dramatically increased 
at time t ~ 11 (from the line length i ~ 32 to t ~ 50) due 
to the rarefaction pulse being transformed into a small 
vortex ring during a short time interval, 5t ~ 5. After 
that the ring loses its energy to Kelvin wave radiation and 
stabilises to axisymmetric vortex ring of slightly smaller 
radius at t ~ 100, so that the total loss of the vortex 
line is about 3%. The length of the vortex line does not 



change for t > 100. We repeated the calculations for 
the rarefaction pulse from the upper branch of the JR 
dispersion curve and found that the ring stabilises after 
losing 4% of its length. These results may seem to con- 
tradict the findings of 6] , where the length of the vortex 
line continues to decay after the collision. The calcula- 
tions of were done in computational box with periodic 
boundaries, so the sound waves generated during the re- 
connections of a large ring with a small ring continue to 
excite Kelvin waves on the vortex line with a continuous 
loss of the vortex line to sound radiation. The rate of 
the decay of the vortex line in this case should depend 
on the size of the computational box. In our calculations 
the faces of the computational box were open to allow 
sound waves to escape pj], which allowed us to compute 
the effect of a single collision with the finite amplitude 
sound wave. 

For very small vortex rings (R ~ 1) the collision with 
a rarefaction pulse causes all the energy and momentum 
to be converted into sound waves with a complete loss of 
the vortex line. 



V. DECAY OF A VORTEX TANGLE 

Next, we estimate the effect that the interactions with 
rarefaction pulses have on a vortex tangle. This time 
we perform the simulations similar to those of 0, IT3 |. 
in which the vortex tangle was created by four collid- 
ing rings in a periodic box. But in our simulations we 
regulate the initial amount of sound waves in our com- 
putational box by introducing them as randomly oriented 
rarefaction waves in addition to the vortex rings. 

For simplicity we shall assume that all sound waves 
initially consist of three types of rarefaction pulses mov- 
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FIG. 7: (colour online) The time snapshots of the density 
isosurfaces of the condensate at 2 = 0.3 for the collision of 
the vortex ring of radius 5.1 and the rarefaction pulse moving 
with the velocity U — 0.63. The density contour plots at 
the cross-section for z — are given below each isosurface. 
Black solid and dashed lines show zeros of real and imaginary 
parts of ip correspondingly, therefore, their intersection shows 
the position of topological zeros. Both low and high density 
regions are shown in darker shades to emphasise intermediate 
density regions. Only a portion of an actual computational 
box is shown. 
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FIG. 9: Density isosurfaces (1^1 2 = 0.3) showing the snap- 
shots of the time evolution of the initial state consisting of four 
vortex rings of radius R — 30.1 and 100 rarefaction pulses. 
Many small vortex rings are formed during the evolution. 
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FIG. 8: (color online) The vortex line length for the collision 
of four vortex rings of radius R — 30.1. Initially the number of 
rarefaction pulses is N rar e — (black line), N ra re — 100 (red 
or dark grey line) and N rar e = 200 (green or light grey line). 
A decay characterized by %2 = 0.3 is plotted (thick line) to 
show the experimentally determined decay of the vortex line 
length. 
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ing in random directions. The approximations to these 
solitary waves, that were developed in 9], become very 
useful in the construction of our initial conditions. If 
the rarefaction pulse given by (J2J is rotated by angle 
a to the x— axis in the (x,y) plane and by angle (3 
to the y— axis in the (y,z) plane, then its wavefunc- 
tion is given by ijj(x',s') where x' = xcosa — ysma, 
s' = \J (y r cos (3 — zsin/3) 2 + (zcos(3 + y' sin/?) 2 , and 
y' = ycosa + xs'ma. The angles of rotation a and (3 
as well as the center of the rarefaction pulse are chosen 
randomly. 

Our calculations were performed in a periodic box with 
volume, V = 80 3 , which is divided by 16 1 3 grid points 



with a spacing 0.5. The numerical scheme is 4th order 
globally accurate with 4th order Runge-Kutta time in- 
tegration with the time step 0.025. A grid spacing of 
0.25 together with the time step 0.00625 was also used 
to test the accuracy of the numerical method. Initially 
there are four vortex rings of the radius R = 30.1 centered 
at (R, 1, 0), (-R, -1, 0), (1, -R, 0) and (-1, R, 0) that are 
moving towards the center of the box. N rare rarefaction 
pulses are distributed randomly in the internal half of the 
box away from the vortices. The results for the dynamics 
of the vortex line length for three sets of computations 
with N rare = 0, 100 and 200 are plotted in Fig. 8. The 
density isosurfaces plots for N rare = 100 are given in 
Fig. 9. In the case of nonzero number of rarefaction 
pulses there is a rapid initial growth of the length of the 
vortex line due to the evolution of some of the rarefaction 
pulses into small vortex rings as the result of the energy 
transfer among them 9]. Consequently, the length of 
the vortex line shows the balance between vortices being 
created and vortex line being destroyed due to reconnec- 
tions and Kelvin wave radiation. Our simulations show 
that apart from the initial growth of the vortex line in 
the simulations with N rare = 100 initially, the dynamics 
for N rare = 100 and N rare = is very similar for in the 
time interval [400,2600]. On the other hand, the case 
with N rare = 200 is quite different. During the time in- 
terval [500, 1300], the decay of the vortex tangle is much 
faster than for two other runs. Whereas after t ~ 1800 
it is the run with N rare = 200 that shows a slower de- 
cay. We speculate that this can be best explained by 
the following. During the moderate times, the vortex 
line loss due to the excitations of the vortex line created 
by rarefaction pulses depends on how many rarefaction 
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pulses survived as their energy is gradually being con- 
verted to high-frequency waves. The more of rarefaction 
pulses were present in the initial state, the longer at least 
some of them are present in the system, therefore, the 
state with N rare = 200 shows a faster decay during the 
moderate times. At larger times, all of the rarefaction 
pulses are destroyed by the collisions and vortices live on 
the background of high-frequency waves that limit their 
mobility and reduce the number of reconnections which 
slows down the decay of the vortex line length, with the 
total energy carried by these high-frequency waves being 
higher for the simulations with N rare = 200 initially. 

For even larger times the situation is reversed one 
more time as reconnections are now infrequent and the 
main decay mechanism is through Kelvin wave radiation. 
Therefore, the rate of decay of vortex line is propotional 
to the amount of high-frequency waves (slowest decay for 
N r are — 0, fastest decay for N rare = 200.) 

For Kelvin- wave cascade, where energy is transfered to 
a much shorter wavelengths with a cutoff below a criti- 
cal wavelength, the vortex line density L = £/V can be 
described by the Vinen equation 



dL k 2 
~dt = -2^ X2L ' 



(8) 



where n = 2tt in our dimensionless units and \2 is a 
dimensionless coefficient. More accurately \2 is a weak 
(logarithmic) function of L and of other parameters such 
as the Kelvin cutoff and temperature. The logarithmic 
dependence on L can be easily obtained in the context 
of the local induction approximation. The presence of 
a large amount of rarefaction pulses changes this weak 
dependence on L and the decay of the vortex line ap- 
proaches an exponential decay instead: 

dL k 

In Fig. 10 we plotted the curve (the solution of (JSJ) 



£ = £ exp(-X2^) 



(10) 



for X2 = 0.00035 and £q = 760 to illustrate the exponen- 
tial decay of the vortex line length. 

Finally, we note that the presence of rarefaction waves 
in a condensate with regular vortex lattices has destabil- 
ising effect and facilitates reconnections. If vortices are 
antiparallel, the Kelvin wave created on the filaments as 
the result of interaction with rarefaction waves facilitates 
the growth of the Crow instability that leads to vortex 
reconnections that destroy a significant amount of vortex 
line (as we estimated in the minimum fractional line 
loss of the pair of antiparallel vortices that are initially 
the distance h apart and perturbe d by a wave of the 
wavenumber k <C 1 is about 1 — y/kh/2). The interac- 
tions of the rarefaction waves with the lattice of parallel 
vortices (created, for instance, in rotating condensates) 
also destabilises the array and creates a tangle as Fig. 11 
illustrates. 



FIG. 10: (color online) The vortex line length for the collision 
of four vortex rings of radius R — 30.1. Initially the number 
of rarefaction pulses is N ra re = 200 (green or light grey line). 
The black solid line gives a decay of the vortex line given by 
JTU| for £ = 760 and X 2 = 0.00035. 
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FIG. 11: Density isosurfaces for \ip\ 2 — 0.3 showing two snap- 
shots of the time evolution of the array of 10 initially straight 
parallel vortices that interacted with 10 rarefaction pulses. 
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VI. CONCLUSIONS 

In summary, we established several mechanisms for 
creation of the rarefaction pulses in condensates. We 
have studied the effects of their interactions with straight 
line vortices, vortex rings and tangles of vortices. We 
showed that there are two competing mechanisms of 
these interactions. Firstly, the interactions of waves at 
close distances lead to the increase of the vortex line as 
the rarefaction pulses in the regions of lower density (in 
the vicinity of the vortex core or each other) may acquire 
circulation and become vortex rings. Secondly, rarefac- 
tion waves excite Kelvin waves on the vortex filaments 
causing a loss of the vortex line due to sound emission. 
At high vortex line densities our simulations suggest that 
the Kelvin wave radiation is enhanced by the presence of 
rarefaction pulses and this accounts for a dramatic in- 
crease in the rate of the decay of the vortex line length. 
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